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Abstract 

Chemical and thermal freeze-out of the hadronic fireball formed in symmetric 
collisions of light, intermediate-mass, and heavy nuclei at beam energies be- 
tween 0.8A GeV and 2. OA GeV are discussed in terms of an equilibrated, 
isospin-symmetric ideal hadron gas with grand-canonical baryon-number 
conservation. For each collision system the baryochemical potential and 
the chemical freeze-out temperature T c are deduced from the inclusive 7r° and 
r\ yields which are augmented by interpolated data on deuteron production. 
With increasing beam energy fiB drops from 800 MeV to 650 MeV, while T c 
rises from 55 MeV to 90 MeV. For given beam energy fig grows with system 
size, whereas T c remains constant. The centrality dependence of the freeze- 
out parameters is weak as exemplified by the system Au + Au at 0.8A GeV. 
For the highest beam energies the fraction of nucleons excited to resonance 
states reaches freeze-out values of nearly 15%, suggesting resonance densities 
close to normal nuclear density at maximum compression. In contrast to the 
particle yields, which convey the status at chemical freeze-out, the shapes 
of the related transverse-mass spectra do reflect thermal freeze-out. The ob- 
served thermal freeze-out temperatures Tth are equal to or slightly lower than 
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T c , indicative of nearly simultaneous chemical and thermal freeze-out. 
PACS numbers: 24.10.Pa, 25.75.-q, 25.75.Dw 
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I. INTRODUCTION 



Relativistic nucleus-nucleus collisions offer the unique possibility to study nuclear matter 
under the influence of high temperature and pressure. At incident energies around 1A GeV 
the formation of a hot and dense reaction zone, dubbed the fireball, has been verified exper- 
imentally M. According to model calculations, the compression phase lasts for time spans 
of 10 - 15 fm/c and reaches densities of 2 - 3 times the nuclear ground-state density [|. 
Simultaneously, temperatures up to ^100 MeV may be achieved and a substantial fraction 
of the nucleons participating in the collision is excited to heavier, short-lived resonance 
states which decay predominantly via meson emission f|-^|. Thus the fireball produced in 
the energy regime of the heavy-ion synchrotron SIS at GSI Darmstadt comprises nucleons, 
resonances, and mesons. 

It is an interesting question to what extent this hadronic system can be described in terms 
of chemical and thermal equilibrium. In the present paper we address this issue on the basis 
of inclusive as well as centrality-selected data on 7r° and 77-meson production in collisions 
of nuclei with equal mass number A. Chapter 2 introduces the relevant concepts and gives 
a brief account of previous studies of the hadronic fireball at SIS energies. The existing 
systematics of 7r° and rj production is reviewed in chapter 3. There we demonstrate that 
thermal concepts do provide a useful description of the fireball at SIS energies. Our model 
of an ideal hadron gas in chemical equilibrium is presented in chapter 4. Chapter 5 discusses 
the results of the thermal analysis. The chemical composition of the fireball at freeze- 
out is determined, and the consistency of chemical and thermal freeze-out temperatures is 
addressed. Finally, the freeze-out conditions of hadronic matter as derived for the SIS energy 
regime are compared to results obtained from similar analyses of particle ratios measured 
at significantly higher incident energies at the AGS (Brookhaven National Laboratory) and 
the SPS (CERN). 
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II. THERMAL CONCEPTS IN NUCLEUS NUCLEUS COLLISIONS AT SIS 

ENERGIES 



In the initial phase of nucleus-nucleus collisions with incident kinetic energies E beam 
around 1A GeV a system of interacting nucleons, resonances, and mesons is created. The 
size of the hadronic fireball depends on the masses of target and projectile nucleus and, 
in addition, on the centrality of the collision. Within the participant-spectator model, 
the geometrical overlap of the two nuclei determines the number of nucleons A part which 
are directly involved. The quantity B = A part thus defines the number of baryons in the 
system. For symmetric collisions and complete stopping the energy available per baryon 
is \fsj1 — mjv, with s given by 4m 2 N + Im^E^am and denoting the nucleon mass. 
The initial conditions of the fireball are therefore fixed. In the subsequent time evolution 
the available energy is transformed into the excitation of thermal and collective degrees of 
freedom. The energy turns into heat and provides the mass stored in resonance states and 
mesons at chemical freeze-out. In addition, the energy builds up compression and produces 
the flow of the expanding matter. 

In a simplified picture hadrons cannot escape from the fireball during the high-density 
phase of the collision. Nucleons, resonances, and mesons are trapped in a cyclic process 
of generation, absorption, and re-emission, exemplified for nucleons N, n mesons, and A 
resonances by NN # iVA # NNir. Within this approach hadrons are released only with 
the onset of the expansion phase, when mesons and baryons decouple due to the decreasing 
matter density. Moments in the expansion process when certain degrees of freedom of the 
system no longer participate in the interaction provide landmarks in the time evolution. 
One has to distinguish between chemical and thermal freeze-out which - in the limit of 
sudden freeze-out - correspond to those moments in time when the relative abundances 
of the particle species or their momentum distributions stop to change. While only inelas- 
tic collisions involving the short-range nuclear force can alter the relative particle yields, 
the momentum distributions of the particles are governed by the larger total interaction 
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cross sections. Consequently, thermal freeze-out does not occur before chemical freeze-out. 
Through the frequent scattering processes of the constituents the system maintains chemical 
and eventually thermal equilibrium. 

It has been shown that at the AGS (E beam < 13.7 A GeV) and the SPS (200A GeV) a large 
number of hadronic observables, including strange-particle yields, is in quite good agreement 
with such an equilibrium scenario j9Hl2j. At the much lower SIS energies, however, only 
a very limited variety of hadron species is produced with significant yields. Therefore, in 
contrast to the situation at AGS and SPS energies, the number of observables which may 
reflect the chemical or thermal freeze-out is comparable to the number of free parameters in 
common thermal models. Nevertheless, several analyses have been performed to check for 
consistency between data and thermal-model predictions also at SIS energies. Midrapidity 
transverse-momentum spectra of charged pions, protons, and deuterons measured in central 
Ni + Ni collisions at 1.06 A 1.45 A and 1.93A GeV were observed to be consistent with 



thermal equilibrium if, in addition, collective radial flow was taken into account [ 13"] . For 
the same reactions also chemical equilibrium has been claimed with close agreement between 
the chemical freeze-out temperatures and the temperatures derived from the particle spectra. 
Recently, Cleymans et al. have extended the thermal analysis at SIS energies to comprise 
also the strange mesons K + and K~. In their systematic study of central collisions of Ni + Ni 



and Au + Au based on a model with canonical strangeness conservation |14],[L5] the yields 
of protons, deuterons, charged pions, and K mesons were found to agree with chemical 
equilibrium at freeze-out. Failure of their analysis to also accomodate the r] meson within 



this common freeze-out picture [15] may be attributed to the use of extrapolated rj yields. 



III. EVIDENCE FOR THERMAL BEHAVIOR FROM MESON YIELDS AND 

SPECTRA 

Pions and 7] mesons are the most abundantly produced mesons at SIS energies. While 
the pionic degree of freedom is also covered by the spectroscopy of the charged members 
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of the isospin multiplet, it is only through 7-ray spectroscopy that the 77 meson becomes 
observable in nucleus-nucleus collisions. The two neutral mesons tt° and 77 can both be 
identified by a two-photon invariant-mass analysis of coincident photon pairs. Using the 
Two-Arm Photon Spectrometer, the TAPS collaboration performed a series of systematic 
meson-production experiments covering incident energies from 0.2A GeV to 2.0A GeV. 
These measurements have established an extended data base for tt° and 77 production in 
light (C + C), intermediate-mass (Ar, Ca + Ca), and heavy symmetric systems (Ni + Ni, 
Kr + Zr, Au + Au) |J-|23 • 



The primary information for the present discussion is provided by the inclusive meson 
yields. They have been determined within narrow rapidity intervals around midrapidity. For 
a consistent treatment the data of the various systems are extrapolated to the full solid angle 
assuming an isotropic source at midrapidity. Fig. [1] shows the resulting average inclusive 
meson multiplicities (M), normalized to the average number of participating baryons (B), 
as a function of the energy available per baryon divided by the corresponding meson mass. 
Except for the system Au + Au at 0.8 A GeV JTB| and the systems Au + Au and Kr + Zr 
at 1.0 A GeV (see table |), where (B) has been determined experimentally, we calculate (B) 
from the geometrical overlap of two colliding sharp spheres which gives (B) = A/2 in case 
of collision partners with equal mass number A. A steep rise of the normalized meson ratios 
with increasing energy is observed. In addition, a weaker dependence on the size of the colli- 
sion system is visible with the clear tendency towards smaller inclusive yields in the heavier 
systems. The same trends have also been observed in charged-pion production experiments 
|24] , |2"5|| . In first approximation all data points fall onto a smooth curve indicating that to 



a large extent the meson-production probability is determined by the energy available per 
baryon. This is quite remarkable because 7r° and rj production proceed through different 
baryon resonances: 7r° mesons mainly come from A(1232)-resonance decays, while the heav- 
ier 77 mesons essentially originate from the iV(1535) resonance which, at SIS energies, is the 
only significantly populated baryon resonance with a large decay width into 77 mesons. The 
fact that the basic production mechanism is no longer apparent in the observed meson yields 



may be interpreted as a first indication for meson emission from an equilibrated source. 

The second important result concerns the average transverse momenta (p t ) at midrapid- 
ity. For the decay of nominal-mass A(1232) and iV(1535) resonances one has center-of-mass 
momenta of 229 MeV/c for tt° from A(1232) and 182 MeV/c for rj from JV(1535). For reso- 
nances embedded in an isotropic fireball one therefore expects average transverse momenta 
obeying (pt) n ° > (Pt)»?- Fig. |] shows the experimental results as a function of the energy 
available per baryon. Although the uncertainties for the rj data are large compared to those 
of the 7T° data, the rj momenta are clearly higher than the tt° values, opposite to the ex- 
pectation. For both mesons the initial rise of (pt) slows down with increasing available 
energy, indicating possible saturation of the average transverse momentum. We conclude 
that also with respect to the shapes of the meson spectra the naive expectations based on 
the production mechanism are not borne out experimentally. 

As a third piece of evidence for thermal behavior we mention the transverse-mass spectra 
of the mesons. At midrapidity the transverse mass m t of a particle is equivalent to its total 
energy in the center-of-mass system. Following P^ , p7fl , the transverse-mass distribution 
of particles with mass m emitted isotropically from a thermal source is characterized by a 
Boltzmann temperature Tb and, at midrapidity, can be approximately described by 



In Fig. m transverse-mass distributions of ir° and rj mesons are plotted together with fits 
according to the parametrization given in Eq. |]for three different systems at incident energies 
near 2A GeV. Within each reaction the ir° and rj spectra exhibit almost identical inverse- 
slope parameters Tb which in addition do not change significantly with the mass of the 
colliding nuclei. In all three systems the 7r° and rj intensities roughly coincide for m t > m v . 
This indicates that it is the energy required to produce a given transverse mass which 
determines the relative abundance of the meson species near midrapidity. For low m t , 
however, individual differences among the collision systems become apparent. Fig. ^| shows 
a systematic enhancement over the exponential rise extrapolated from the high-mj region 
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(1) 



if one goes from the light C + C to the heavy Ni + Ni system. The same observations, 
namely m t scaling of the 7r° and 77 intensities and low-m t enhancement of the 7r° spectrum 
in heavy systems, have been reported at various energies down to 0.8A GeV and seem 
to be a general feature of heavy-ion collisions in the SIS-energy regime [p^ . p!7| , p!9| - |2T| , |23 



Possible explanations that have been suggested for the low-m t enhancement involve pion 



rescattering through resonance states in the heavier systems and multiple pion decay of 



heavy resonances [p9 



Transverse-mass scaling in the high-m 4 region does not prove that the fireball actually 
has reached chemical equilibrium. The particle spectra of an equilibrated source, however, 
would follow phase space distributions as observed, at least for vanishing flow. 



IV. IDEAL HADRON GAS MODEL 

In our ansatz we assume that at chemical freeze-out the fireball can be described in 
terms of an ideal, equilibrated hadron gas. As constituents we take into account pions and rj 
mesons, as well as nucleons, deuterons, and all non-strange baryon resonances up to a mass 
of 1.8 GeV (see table HD, which corresponds to yfs — m N at 2 A GeV beam energy. 

In the grand-canonical description chosen, hadronic matter is characterized by a baryo- 
chemical potential /i# and a temperature T c . Furthermore, one might consider the isospin 
of the system. In the present analysis, however, we neglect isospin as an additional degree 
of freedom, mainly because our data base of 7r° and rj yields is insensitive to isospin. Within 
the isobar model for instance, the ir° yield amounts to one third of the total pion yield of 
a heavy-ion collision irrespective of the actual isospin of the projectile-target system. The 
other important observable at SIS energies, the r] meson, does not carry isospin. Thus, only 
the third observable in our data base, the relative yield of deuterons and nucleons, depends 
on the isospin. 

In the grand-canonical description of a system of noninteracting fermions and bosons 
the particle-number densities are given by integrals over the particle momentum p 
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where Qi is the spin-isospin degeneracy, the total energy in the local restframe, and Bi 
the baryon number of the particle species i. The form of the denominator in the integrand 
accounts for the different statistics of fermions and bosons. 

Eq. ^| describes an infinitely large system of stable particles and cannot be directly 
applied to the hadronic fireball in a nucleus-nucleus collision. To account for the fact that 
the fireball occupies only a finite volume, a surface correction has to be included. We assume 
a spherical freeze-out volume V c of radius R c and obtain a momentum-dependent correction 
factor f(pR c ) 0,0] 



1 

in the integrand of Eq. 0. The correction leads to a 30-40% reduction of the individual 
particle-number densities as compared to infinite nuclear matter. Particle ratios, however, 
are hardly affected since here the corrections nearly cancel. The freeze-out radius of the 
system is fixed by baryon number conservation 

/ PB dV=(B), (4) 

where the baryon density ps comprises contributions from nucleons, deuterons, and reso- 
nances 

Pb = Pn + 2pd + ^2 Pr- (5) 

R 

Given the feedback via Eq. |^, R c has to be adjusted iteratively. Starting from an ad-hoc 
initial value, R c is varied in steps of 0.5 fm in order to fulfill Eq. |] within the uncertainties 
of the experimental particle ratios. 

The second modification of Eq. |2] is related to the fact that the baryon resonances are 
unstable particles which do not have a fixed mass but exhibit a broad mass distribution. 
In contrast to thermal-model analyses of AGS and SPS data, the actual mass distributions 
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cannot be neglected in the 1A GeV energy regime since here the energy available in the 
nucleon-nucleon system is comparable to the excitation energy of the baryon resonances. 
This means that the A(1232) resonance, the lowest of the baryon resonances, is by far the 
most abundantly populated resonance and that all resonances are populated predominantly 
in the low-mass tail of their mass distributions. To account for the off-shell behavior of the 
resonances, the integrand in Eq. ^|is folded with normalized mass distributions Ai(m). These 
are trivial delta functions for the stable nucleon and deuteron and for the electromagnetically 
decaying ir° and 77, while the scheme to parametrize the resonance mass distributions has 



been adopted from [29 



We describe a given resonance mass distribution by a Lorentzian of the form 

m 2 T n(m) 



A R( m ) °c T-r, ""V^ 9/ \ , (6) 



(m 2 — ra R ) 2 + m 2 r| 

with 



R\ m > 



00 

A R (m) dm = 1 . (7) 

The characteristic parameters are the nominal resonance mass m R and the mass-dependent 
total decay width T R (m). As decay modes we take I71", 2tt, and r\ decay to the nucleon 
ground state into account. The total decay width at mass m thus is given by 

r R (m) = r l7r (m) + r 2vr (m) + r„(m) . (8) 

The 27r decay is described in terms of a one-step process where the resonance decays into a 
nucleon and an object with angular momentum / = and twice the pion mass = 2m 7r . 
The hypothetical di-pion subsequently disintegrates into two pions. With this simplification 
the 27r decay width can be written in complete analogy to the In and rj widths. For the 
mass-dependent partial decay widths r l7rj 27r,jj(m) we thus obtain 

, J i\ 2l+1 (q 2 R + t 2 \ l+1 

T ln ,2n, v ( m ) = { V<W W + 52 J (9) 

otherwise 
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The partial decay widths at m = m# are listed in Tab. ||. The quantity I is the angular 
momentum of the emitted meson or di-pion, and q is the momentum of the particle in the 
restframe of the decaying resonance, with denoting the special case of m = hir. The 
cutoff parameter 5 is given by 



2 r 2 [rriji) 

5 2 = (m R - m N - m-n^^f + - j - • (10) 

Only in the case of the A(1232) resonance we deviate from this prescription, mainly for the 
sake of consistency. Use of the well-established parameters Tr = 110 MeV and 5 = 300 
MeV/c determined for this resonance by Koch et al. in a non-relativistic analysis (see P3| ) 
requires an additional term hir/hi on the right-hand side of Eq. [| 

Furthermore, one could consider an excluded-volume correction to take into account the 
hadron-hadron hard-core repulsion. This means that one would consider a real rather than 
an ideal hadron gas. An excluded-volume correction, however, does not play a significant 



role for particle-yield ratios as demonstrated in |12|] . Another reason to omit the excluded- 
volume correction is the fact that chemical freeze-out seems to occur at baryon densities 
well below the nuclear ground-state density. The system therefore is rather dilute which 
supports the assumption of an ideal gas of pointlike non-interacting particles. 

With the modifications introduced by Eqs. || and || the particle-number densities of Eq. |2| 
are given by 



Pi = ^2 I d PP 2 f(P R c) 
2tt Jo 

Ai(m) 



dm ■ 



n 



± 1 



(11) 



exp (y/m 2 + p' 2 - HbB^j /T c 
Equation ^ employs a mass- dependent width in the resonance-mass distribution. This fea- 
ture of the model makes it possible to take the actual resonance masses into account if one 
determines the contributions to the asymptotically observed 7r° and rj intensities which orig- 
inate from resonance decay. The effective branching ratios may deviate significantly from 
those at the nominal mass hir. For the latent meson densities p]i ,27r,v represented by a given 
resonance R one obtains 
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l7r,27r,»7 
Pr 



W 



9r 

2tt 2 J u 
dm ■ 



dpp 2 f(pR c ) 

^(i7T,2n^)(m)A R (m) 



(12) 



'o T R (m)exp U/rri 2 + p 2 — ^b) JT C +1 

with w — 1 in the case of one-pion or ^-meson emission and with w = 2 for the two-pion 
decay. The resonances present at freeze-out also contribute to the asymptotically observed 
abundance of nucleons. Furthermore, the i] meson (lifetime cr = 1.7 x 10 5 fm) gives rise 
to additional 7r° intensity via its 37r° and 7r + 7r~7r° decays with branching ratios of 32% and 
23%, respectively. 



V. THERMAL MODEL ANALYSIS 

The values for the two parameters of the hadron-gas model - the baryochemical potential 
fiB and the temperature T c - can be derived from any two relative abundances of the con- 
stituents. Furthermore, the baryon content of the system has to be known for an assessment 
of the geometrical size of the fireball. The n° and ^-meson yields are particularly suited for 
such an analysis, as these particles unambiguously arise from the fireball. Hence, we base 
our analysis on the ratios (M n o)/ (B) and (M v ) / (M n o) as compiled in Tabs. | and |TTT| . For 



the sake of redundancy in the analysis we consider the ratio of deuterons to nucleons as 
additional input information. The knowledge of (M^)/(M-^), however, is limited and espe- 
cially its centrality dependence is difficult to determine experimentally since in non-central 
collisions it is not trivial to separate deuterons originating from the fireball from those which 
are emitted by target- or projectile-like spectator remnants. 

So far, only a few investigations have addressed deuteron yields in nucleus-nucleus colli- 
sions at SIS energies. The inclusive values for (Md) / (Mn) given in Tab. | are based on yield 
ratios of deuterons and protons measured in central collisions of Ni + Ni at 1.06, 1.45, and 



1.93A GeV [O] and of Au + Au at 1.CL4 GeV |3~5|, and on inclusive results obtained for 



various combinations of light and intermediate-mass nuclei at 0.8A GeV ||36|| . ^From these 



data we parametrized the dependence of (M d )/(M N ) on the incident energy and on the 
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system size. For given system size a linear decrease of (M d ) / (Mn) with increasing available 
energy provides a fair description of the existing data. For fixed energy the ratio (Md) / (Mn) 
increases with system size until saturation sets in for heavy systems. The limited number of 
data points and the uncertainties in the interpolation procedure translate into quite large 
errors for the ratios (M d )/(M N ) (see Tab. |). No attempt has been made to establish an 
impact-parameter dependence. 

In the thermal analysis the three particle ratios are related to the particle-number den- 
sities calculated within the thermal model by the following equations: 



(B) 



exp 



+ 1.2(p, + ^ (1535) )} , (13) 



(M^) \ P v + P^(l535) 



R 



' (M d ) \ p d 

, ( M N) / exp PN + Y.PR 

R 



(14) 



(15) 



Each equation defines a band in the /iB,T c -diagram. The widths of these bands are given 
by the uncertainty of the input values, and for each individual collision system which is 
characterized by its triplet of particle ratios one obtains the chemical potential and the tem- 
perature as coordinates of the common intersection point of the three bands. Provided the 
experimental uncertainties correspond to Gaussian errors one has a probability of 20% that 
the true values of the three particle ratios simultaneously fall into the la uncertainty ellip- 



soid, while this probability rises to 75% for 2a uncertainty [|37|. For a consistent treatment 
of the full data sample it is therefore necessary to allow for uncertainties in excess of two 
standard deviations in the input ratios. 

The system size also has to be considered as input information. Through baryon conser- 
vation expressed in the grand-canonical form of Eq. [| the baryon number determines the 
freeze-out radius R c . For typical values of R c around 5 fm we find that changes within ±1 fm 
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of the final solution hardly affect the numerical values for ps and T c , but substantially alter 
the freeze-out density ps due to its dependence. Larger variations up to ±2 fm may 
still lead to solutions on the 2cr-level for the particle ratios with changes in ps and T c within 
±5%, but ps becomes increasingly incompatible with Eq. [|. 

The connection between a given particle ratio and the freeze-out parameters fig and 
T c is illustrated in Fig. ^. A global value of 5 fm for the radius parameter R c has been 
chosen to generate the diagram. Curves of constant particle ratio in the pB,T c -p\ane are 
referred to as freeze-out lines. Since R c is fixed, (B) varies along the freeze-out lines. The 
intersection of any two freeze-out lines from different particle ratios thus only describes a 
physical solution if also the corresponding values of (B) overlap. The third freeze-out line 
then has to intersect within the experimental uncertainty and for the same value of (B). In 
Fig. [| global values for the triplet of particle ratios have been chosen to demonstrate the 
variation of the freeze-out lines with the radius parameter. 

A. Centrality Dependence of the Freeze out Parameters 

The inclusive particle ratios of the existing data base represent averages over the impact 
parameter of the collision. Prior to the analysis of the inclusive data we therefore address 
the centrality dependence of the freeze-out parameters. 

The ratios (M n o) / (B) and (M v )/(M T o) have been measured as a function of the cen- 
trality of the collision in Ar + Ca and Au + Au |16|, both at 0.8 A GeV. While in 
the intermediate-mass Ar + Ca system the meson multiplicities per participant nucleon 
essentially do not depend on the impact parameter, a clear effect is observed in the heavy 
system Au + Au. ^From peripheral ((B) = 40 ± 10) to central collisions ((B) = 345 ± 25), 
the neutral-pion multiplicity per participant nucleon (M n o)/(B) increases by a factor of 
about 1.6, while simultaneously the ^-meson multiplicity per participant nucleon (M V )/(B) 
increases by a factor of 3.6 (see Tab. pTT] ). The enhancement of the 77-meson production in 
central collisions has been attributed to secondary reactions involving resonances as inter- 
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mediate energy storage ||16|| . This interpretation is particularly relevant for the large system 
size and the subthreshold incident energy of the measurement and indicates that dynamical 
effects may still be visible in these Au + Au collisions. 

Fig. H illustrates the thermal analysis of Au + Au for peripheral, semicentral, and central 
collisions. Shown are the three bands in the /iB,T c -plane corresponding to the three particle 
ratios for la and 2a uncertainties. It is important to understand the correlations between the 
particle ratios and the parameters of the thermal model. One observes that in the vicinity 
of the common intersection point of the bands the ratio (M n o)/(B) is equally sensitive to 
fiB and T c . Larger values of (M 7V o)/(B) mean larger temperatures and smaller chemical 
potentials of the hadron gas. This reflects the fact that with increasing /ig the baryon 
density grows while the density of free pions remains unaffected. The ratio (M v )/(M n o), 
on the other hand, is practically only sensitive to T c . Larger values of (M v )/(M n o) require 
larger values for T c , while hb is uncritical. This is due to the fact that producing an rj meson 
or exciting an iV(1535) resonance requires more energy and thus a higher temperature than 
producing a pion or exciting a A (1232) resonance because of the pronounced mass difference 
between both mesons or both resonances, respectively. The third ratio, (Md)/ (My), derives 
its sensitivity to \lb from the difference in baryon number of the two particles, while the 
temperature dependence is governed by the mass difference between the particles. With 
increasing baryochemical potential the temperature of the system has to drop in order to 
maintain a given (MA / (Mpj) ratio. It should be noted, however, that the sensitivity of 
the freeze-out parameters with respect to a variation of (M^)/(M^) is small. Modifying 
the ratio by as much as ±20% changes the freeze-out parameters by not more than 5%. 
Nevertheless, (Mj)/(Mn) provides a useful additional constraint in those cases where the 
neutral-meson data alone would leave room for large variations in \ib and T c . 

Fig. |] shows that in 0.8 A GeV Au + Au the experimental particle ratios do define unique 
sets of hb and T c for each of the three impact parameter selections. The numerical values 
as determined by x 2 minimization are summarized in Tab. |TJ. The uncertainties quoted 
represent la standard deviations and reflect the size of the error ellipse of the hb,T c — pair 
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at X 2 — Xmin + 1- Within these uncertainties the baryochemical potential and the chemical 
freeze-out temperature obtained in the most peripheral and central collisions, respectively, 
agree with the semicentral values of ps = 812 ± 5 MeV and T c = 52 ± 2 MeV. The ensuing 
baryon densities ps relative to the nuclear ground-state density po are also given in Tab. [[II]. 

In our ansatz we are able to describe the central Au + Au collisions at 0.8A GeV with 
one common set of freeze-out parameters, see Fig. 0, while Cleymans et al. |15| emphasize 



that they cannot reproduce the ratio (M v )/{M n o) within their hadron-gas model for the 
same system at 1.0 A GeV. The present analysis is based on the directly measured ratio 
(M r? )/(M n .o) for experimentally determined values of (B), while the authors of ref. [ll| have 
to overcome the difficulty that at 1.0 A GeV the ratio (M v )/(M n o) has not been measured 
for truely central collisions. They extrapolate (M v )/(M n o) measured in centrality-biased 
Au + Au collisions at 1.0 A GeV beam energy [|17[] (c.f. Tab. |) to fully central collisions 
using the centrality dependence of that ratio measured at 0.8 A GeV [plj. It is known |T6| , |2"4 



however, that the dependence of the meson multiplicity on the number of baryons in the 
fireball is governed by the transverse-mass excess (m t ) — (y^ — 2mAr), giving rise to a steeper 
increase of the r) multiplicity with centrality in the more subthreshold Au + Au collisions at 
0.8A GeV. With their approach Cleymans et al. thus overestimate the (M v ) / (M n o) ratio since 
at 1.CL4 GeV the production of rj mesons is considerably less subthreshold than at 0.8A GeV. 
In addition, Cleymans et al. were not aware of the fact that the 1.0 A GeV measurement is 
already centrality-biased, the value of (1.4 ± 0.6)% corresponding to (B) = 164 ± 20 (c.f. 
Tab. |). Judging from the full variation of 40% observed for the centrality dependence of 
(M 7) )/(M 7r o) at 1.0A GeV (see ref. |T?j) we obtain (M^)/(M 7r o) = (1.8 ± 0.7)% for central 
collisions. Thus, the lower boundary of the 2a uncertainty band for a revised 77/71"° freeze- 
out line may reach T c = 55 MeV at ps — 800 MeV (see Figs. |] and ||), well within the 2a 



ellipsoid of the intersecting remaining freeze-out lines in the analysis of Cleymans et al. [Jlo 

In summary, no significant centrality dependence of the freeze-out parameters is found 
for the heavy system Au + Au at 0.8A GeV. For the lighter systems and for beam energies 
which are less subthreshold for 77 production a possible dependence on the impact parameter 
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therefore can be neglected at the present level of accuracy. 



B. Results for Inclusive Collisions 

Fig. |8] shows the results of the thermal analysis for the lightest, an intermediate-mass, 
and the heaviest system, both at the lowest and highest beam energy studied. Presented 
in the /i£,T c -plane are the bands defined by the input values (M^o) / (B), (M v ) / (M n o), and 
(Md) / (Mn). In all cases the three bands overlap and, therefore, define the freeze-out param- 
eters fiB and T c in an unambiguous way. Even collision systems with as few participating 
nucleons as C + C seem to comply with the model-assumption of chemical equilibrium. The 
results for the full set of inclusive measurements are summarized in Tab. [TV]. The analysis 
reveals a systematic reduction of the baryochemical potential from 800 to 650 MeV with 
increasing beam energy, which is accompanied by an increase of the freeze-out temperature 
from 55 to 90 MeV. For given beam energy the chemical potential /xg grows with increas- 
ing system size while the freeze-out temperature T c stays almost constant. These findings 
agree with results quoted from similar analyses [|T^-[T3| , although we observe a general trend 



towards smaller baryochemical potentials and higher temperatures as compared to | P^UI5" 

The freeze-out radius R c is fixed by the average number of baryons through Eq. |] and 
enters the model analysis via the surface correction term. R c is larger for the heavy target- 
projectile combinations, but increases less than the trivial (B) 1 ^ 3 law, and it decreases with 
increasing beam energy. The corresponding baryon densities at chemical freeze-out are in 
general smaller than about half the nuclear ground-state density which is in good agreement 
with results from other analyses flO|,|TT|,|i"3| . 



C. Chemical Composition of the Fireball 

The quantities fiB and T c are the two important free parameters of our model. They 
essentially determine the properties of the fireball. As an example, Fig. |9| shows the chemical 
composition of the baryon sector at chemical freeze-out as a function of the energy available 
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in the nucleon-nucleon system. The system size is a further parameter. Light, intermediate- 
mass, and heavy systems are therefore treated separately in Fig. [5[ 

Immediately apparent is the pronounced energy dependence of the chemical composition, 
while for any given incident energy the chemical composition is nearly independent on the 
size of the system. The fraction of baryons excited to resonance states grows from 2-3% 
at 0.8A GeV beam energy to about 15% around 2.0A GeV. The A(1232) resonance is 
populated most abundantly. The ratio of heavier resonances to the A(1232) resonance, 
however, increases from about 7% at 0.8A GeV to nearly 20% around 2.0A GeV, indicative 
of a more equal population of the resonance spectrum at higher incident energies. As 
illustrated by the straight lines in Fig. |9|, the energy dependence of the individual relative 
baryon populations is in reasonable agreement with exponential behavior. The slopes of 
these exponentials are nearly equal for A (1232) and iV(1535) and steeper for the summed 
contributions of the remaining A and N resonances. 

The baryon composition of the chemical freeze-out state being established it is interesting 
to extrapolate back to the high-density phase of the collision. Around 2A GeV beam 
energy microscopic model calculations quite consistently give a maximum baryon density 
of 

Pmax ~ 2.5 po an d predict a ratio of 0.3 to 0.4 for the ratio of the number of A(1232) 
resonances at freeze-out to the corresponding number at maximum compression p9| , |38| . 
Using this information together with an allowance for contributions from higher resonances, 
one obtains a value of 0.4 p max ~ po for the maximum resonance density in the collisions at 
1.9 A and 2. OA GeV, see Fig. ||. Although the density of baryon resonances only amounts to 
about 40% of the total baryon density, resonance-resonance interactions might take place. 
Hadronic matter in that state has been referred to as resonance matter ||||. 

Concerning the meson sector of the freeze-out state, an important observation is that 
a sizeable fraction of the mesons are present as free mesons in chemical equilibrium with 
the baryons. This is immediately apparent, if one compares the population of A(1232) and 
A^(1535) resonances with the ir° and r] multiplicities observed asymptotically (see Tab. |). 
For a quantitative discussion, we plot in Fig. |X0| the fraction of free mesons and the fraction 
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of mesons bound in resonances at freeze-out as a function of the energy available in the 
nucleon-nucleon system, again differentiating between light, intermediate-mass, and heavy 
systems. The relative contributions to the final-state pion yields reveal a moderate energy 
dependence, the fraction of free pions decreasing from about 65% at 0.8A GeV beam energy 
to about 55% at 2.0A GeV. Pions from the A(1232) resonance behave in a complementary 
way, as expected, with contributions of about 30% at 0.8A GeV and about 40% at 2.0A GeV. 
The remaining intensity can be attributed to heavier resonances which, according to their 
higher mass, give rise to a steeper energy dependence than exhibited by pions from A(1232)- 
resonance decays. For rj mesons resonance decay after chemical freeze-out is unimportant. 
About 85-90% of the asymptotically observed mesons are already free at chemical freeze- 
out. The energy dependence of the iV(1535) resonance contribution is comparable to that 
of the pion contribution from A(1232) decays (see also Fig. |9[). For pions as well as for rj 
mesons a trend towards larger contributions from resonance decays is visible in the heavier 
systems. 



D. Time Order of Chemical and Thermal Freeze out 

In a sudden freeze-out scenario the temperature of an equilibrated hadron gas at thermal 
freeze-out is characterized by the momentum spectra of the emitted particles. The experi- 
mental spectra, however, may be modified by resonance decays and, in particular in heavy 
systems, by the rescattering of particles off spectator material. In addition, energy trans- 
ferred into collective flow reduces the temperature of the system. Thus, in order to extract 
the freeze-out temperature T th from the spectral shapes, further assumptions have to be 
made leading to model-dependent results. Nevertheless, the meson spectra of the present 
data base do provide a valuable consistency check for the chemical freeze-out analysis in 
the sense that the different model temperatures are subject to constraints due to the time 
order of chemical and thermal freeze-out. 

As an example we consider the impact-parameter inclusive transverse-mass spectra of 
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7T and 7] mesons presented in Fig. |3|. The spectra can be described reasonably well by 
exponential distributions, provided one excludes the low transverse masses from the fit. 
For all three collision systems the Boltzmann parameters Tg, which reflect the temperature 
T t h of the fireball at thermal freeze-out in case of vanishing collective flow, agree with 
the corresponding chemical freeze-out temperatures T c . A similar quality of agreement is 
observed for Ar + Ca at 1.5A GeV, see Tab. |V]. At the lower beam energies of 1.0A and 
0.8A GeV the experimental situation for Tb seems unclear at present, except for the heavy 
systems Kr + Zr and Au + Au. In Au + Au at O.&A GeV we have Tb > T c at the la level, 
while at 1.0 A GeV the Boltzmann temperatures in both systems are significantly higher 



than the corresponding chemical freeze-out temperatures, see Tab. [IV. 

The absence of a low-m 4 enhancement in the C + C spectrum of Fig. |3| is remarkable. 
To the extent that rescattering can be assumed negligible, the light C + C system should 
show the expected influence of resonance decay. Instead, a perfectly exponential behavior is 
observed over the full m t range, although pions from resonance decay do comprise 40% of the 
total pion intensity (see Fig. [TD|) . Consequently, the onset of the low-m t enhancement visible 
in the Ca + Ca and Ni + Ni cases of Fig. |3| probably has to be attributed to rescattering in 
spectator material. 

Collective flow affects the meson spectra. The point is, that one has to know the un- 
derlying flow profile in order to analyze a given spectral shape. For the present analysis 



we have chosen the blast model proposed by Siemens and Rasmussen ||39|| . In this model, 



the fireball in thermal equilibrium is assumed to expand isotropically. All particles in the 
fireball share a common temperature T$r and have a common radial-flow velocity fisR- The 
modification of the spectra compared to pure Boltzmann distributions becomes more sig- 
nificant the heavier the considered particle species is. Since pions and rj mesons have low 
masses their spectra are not very sensitive to the radial-flow velocity. Therefore, a fit to the 
mesonic transverse-mass spectra considering both Tsr and (3sr as free parameters would 
determine the freeze-out parameters only with large uncertainties. To avoid this situation, 
we exploit the fact that for the present range of collision systems (3sr is known to have values 
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between 0.25 and 0.35 in accord with the systematics of radial flow velocities measured for 



heavier particles (A=l to 4) in central collisions as quoted in [ 4"0|| . Thus only Tsr is treated 
as a fit parameter while (3sr is kept constant. In the 7r° cases we furthermore restrict the 
fit to rrit > 400 MeV as was done for the Boltzmann fits. Under these conditions the 7r° 
and 7] spectra can be described within the Siemens-Rasmussen model for all beam energies 
and target-projectile combinations. In general, the extracted temperatures Tsr are 10-20% 



smaller than the corresponding Boltzmann temperatures Tr, see Fig. [11]. The actual re- 
duction is correlated with the magnitude of Psr chosen in the fit, with Tsr decreasing for 
increasing (3 sr. 

Thermal freeze-out does not occur before chemical freeze-out and, therefore, the thermal 
freeze-out temperature T t h cannot be larger than T c . The available midrapidity spectra of the 
tt° and r\ mesons do support this conjecture. If one neglects radial flow, the spectral shapes 
observed in the light and intermediate-mass systems provide inverse-slope parameters which 
are in accord with the chemical freeze-out temperatures T c derived from the particle yields. 
Taking radial flow into account results in slightly lower temperatures at thermal freeze-out 
compared to chemical freeze-out, indicative of cooling as the system develops in time. The 
meson spectra observed in the heavy systems Kr + Zr and Au + Au, in particular, do 
require the inclusion of flow in order to achieve consistency (T c > T t h) between thermal and 
chemical freeze-out temperatures. 

E. Chemical Freeze out Curve for Hadronic Matter 



Fig. [12] shows our results within a schematic phase diagram of hadronic matter, together 
with other data points obtained from particle-production experiments at higher energies 
I0| , [TT[| . In contrast to the AGS and SPS results the chemical freeze-out parameters deduced 



at SIS energies are far below the expected phase boundary between hadron gas and quark- 
gluon plasma. For zero chemical potential the critical temperature is constrained by lattice 
QCD calculations. For finite values of \ir the phase boundary can be approximated in a 
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simple model by equating chemical potential and pressure in the hadronic phase and in an 
idealized quark-gluon plasma JTTJ]. It has recently been noticed by Cleymans and Redlich 
that, within their hadron-gas model, the chemical freeze-out curve corresponds to an average 
energy per hadron of 1 GeV pEIfl . The resulting curve is shown in Fig. [12] as a solid line 
which is in close agreement with the deduced freeze-out parameters. 

At SIS energies thermal and chemical freeze-out seem to nearly coincide. In contrast, 
at SPS beam energies the freeze-out parameters clearly indicate that thermal freeze-out 
reflects a later stage of the collision as significantly lower temperatures compared to chemical 



freeze-out are deduced [42 



VI. SUMMARY 

The inclusive neutral-pion and //-meson yields measured in symmetric collisions of light, 
intermediate-mass, and heavy nuclei near midrapidity in the energy range from 0.8A GeV to 
2. OA GeV are consistent with the formation of a hadronic fireball in chemical equilibrium, as 
described by an isospin-symmetric ideal hadron gas. With increasing bombarding energy the 
baryochemical potential \ib decreases from 800 MeV to 650 MeV while simultaneously the 
temperature T c increases from 55 MeV to 90 MeV. Concerning the system-size dependence, 
we find that hb grows with increasing mass of the projectile-target combination while T c 
remains about constant. The centrality dependence of \xb and T c has been investigated in 
the system Au + Au at 0.8A GeV. Here the rj/ir° ratio is expected to be most susceptible 
to the impact parameter, both because of the large mass of the system and because of the 
low incident beam energy. No centrality dependence is observed. 

Apart from the moderate system size dependence, the freeze-out parameters fis and T c 
completely characterize the system. Given \xb and T c one can calculate the hadrochemical 
composition of the fireball. While at 0.8A GeV beam energy only 2-3% of the nucleons are 
excited to resonance states at chemical freeze-out, this fraction increases to about 15% at 
2. OA GeV bombarding energy. According to transport-model calculations, this resonance 
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content at freeze-out implies a resonance density in the high-density phase of the collision 
of about normal nuclear matter density, thus justifying the term resonance matter for the 
specific state of hadronic matter created in heavy-ion collisions at SIS energies. In the 
meson sector of our model about 50% of all asymptotically observed pions are still bound in 
resonance states at chemical freeze-out, while only about 10% of the final-state rj mesons 
are bound at that moment. 

In contrast to the particle yields, which convey the status at chemical freeze-out, the 
shapes of the related transverse-mass spectra do reflect thermal freeze-out. After proper 
allowance for radial flow the slope parameters observed for the midrapidity spectra of the 
neutral mesons correspond to thermal freeze-out temperatures which are equal to or slightly 
lower than the chemical freeze-out temperatures. In contrast to ultrarelativistic collisions 
studied at the SPS, chemical and thermal freeze-out thus seem to occur almost simultane- 
ously at SIS energies as is also observed for heavy-ion collisions at the AGS. 

Particle yields and spectral shapes reflect only two facets of the very complex process 
of a relativistic nucleus-nucleus collision. With these observables alone it is not possible to 
decide whether chemical equilibrium is actually reached during any stage of such a collision. 
In fact, the observation that for the heavy collision systems the agreement between data 
and hadron-gas model is not as good as for the light systems can be seen as one indication 
that physics beyond equilibrium concepts may still be visible in the final state of relativistic 
nucleus-nucleus collisions around 1A GeV. 
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FIG. 1. Average tt° and 77-meson multiplicities per average number of participants in symmetric 
nucleus-nucleus collisions as a function of the energy available per baryon normalized to the mass of 
the respective meson. The data are taken from [p^-p3|. The curve represents a fit to the data and is 
given by the polynomial expression log((M) /(£?)) = —2. 102+ (3. 25— (1.405+0.785 log x) log x) log x, 
where x is the normalized meson-specific energy available. 
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FIG. 2. Average 7r° and 77-meson transverse momenta measured in narrow rapidity intervals 
around midrapidity as a function of the energy available per baryon. The data are taken from the 
same experiments as listed in Fig. [|. The data points at available energies of 182 MeV (0.8,4 GeV 
beam energy), 223 MeV (1.0 A GeV), and 410 MeV (2.0,4 GeV) are slightly shifted in energy with 
respect to their nominal position (center of brackets) to make the error bars of the r] momentum 
visible. With the exception of Ne + Al the error bars for the 7r° momentum are smaller than the 
symbol size. Parallel lines interpolating the 7r° and 77 data, respectively, are drawn to guide the 
eye. 
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FIG. 3. Impact-parameter inclusive transverse-mass spectra of ir° and ry-mesons as observed 



in the systems C + C @ and Ca + Ca §3fl at 2.0 A GeV beam energy and in Ni + Ni at 1.9A GeV 
[19]. The distributions are divided by the square of the transverse mass. In this representation 
midrapidity particles from a thermal source are expected to exhibit a purely exponential spectrum. 
The solid lines represent Boltzmann fits (see Eq. ||) to the 7r° data for mt > 400 MeV. 
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FIG. 4. Freeze-out lines in the /Us , T c -plane for selected values of the particle ra- 
tios 7r°/B, T]/ir°, and d/N. The diagram corresponds to a fixed size of the system with 
the consequence that the baryon density pb(pb,T c ) changes. Representative densities are 
Pb(650 MeV, 70 MeV)=0.053 p , Pb(650 MeV, 100 M e y)=0.469 p , pu(800 MeV, 80 MeV)=1.197 
P0j and /?s(800 MeV, 50 MeV)=0.125 po, where is the nuclear ground-state density po = 
0.168 mi" 3 . 
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FIG. 5. Freeze-out lines in the /xb, T c -plane for selected values of the freeze-out radius R c . The 
diagram corresponds to fixed particle ratios 7T°/B, r]/ir°, and d/N, while R c varies in equidistant 
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FIG. 6. Determination of the chemical freeze-out parameters in the system Au + Au at 
0.8 A GeV beam energy for peripheral ((B) = 40 ± 10, left), semicentral ((B) = 227 ± 20, middle), 
and central collisions ((B) = 345 ± 25, right). For given impact parameter the particle ratios 
7r°/B, rj/ir°, and d/N define three bands, shown in different greyscales, in the T c -plane. The 
solid (hatched) bands reflect la (2a) intervals in the experimental uncertainties. The freeze-out 
parameters are determined by the overlap of the three bands. 
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FIG. 7. Comparison between experimental particle ratios and ratios as calculated from the 
baryochemical potential \ib and the temperature T c at freeze-out. The case shown corresponds to 
central Au + Au collisions at GeV with baryon number (B) = 345 ± 25. R c is the freeze-out 
radius. 
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FIG. 8. Determination of the chemical freeze-out parameters T c and \ib- For given system 
and beam energy the particle ratios vr°/B, t//vt°, and d/N define three bands, distinguished by 
different greyscales, in the hb, T c — plane. The solid (hatched) bands reflect la (2a) intervals in the 
experimental uncertainties. The freeze-out parameters are determined by the overlap of the three 
bands. The individual panels show the solutions obtained for the lightest, an intermediate-mass, 
and the heaviest system studied at the lowest and highest beam energies, respectively. 
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FIG. 9. Baryon content at chemical freeze-out in the hadron-gas model. The relative contri- 
butions to the total baryon number from nucleons and deuterons, A(1232) resonances, 7V(1535) 
resonances, and the sum of all remaining A and N resonances are plotted as a function of the 
energy available per baryon. The three panels show from left to right the results obtained for light, 
intermediate-mass, and heavy systems. The error bars reflect the uncertainty in the freeze-out 
parameters [is and T c (see Tab. |Tv| ). Straight lines are drawn to guide the eye. 
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FIG. 10. Pion (left) and r\ meson sources (right) at chemical freeze-out in the hadron-gas 
model. The mesons are either present as free particles or still bound in resonance states. The 
corresponding relative contributions are shown as a function of the energy available per baryon 
for light (top), intermediate-mass (middle), and heavy systems (bottom). The error bars reflect 
the uncertainty in the freeze-out parameters \xb and T c (see Tab. fV\ ). Straight lines are drawn to 
guide the eye. 
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FIG. 11. Chemical (T c ) and model-dependent thermal (Tg and Tsr) freeze-out temperatures 
as a function of the energy available per baryon. The data points at 0.8A GeV, 1.0 A GeV, and 
2.0^4 GeV beam energy are slightly shifted in energy with respect to their nominal position (center 
of brackets) . 
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FIG. 12. Phase diagram of hadronic matter. The chemical freeze-out temperatures T c are 
shown as a function of the corresponding baryochemical potentials /ig as obtained from yield 
ratios of particles produced in nucleus-nucleus collisions at various incident energies. Results of 
the present work are plotted together with AGS and SPS results 1Q|, The grayscale of the SIS 
data points reflects the system size, i.e. light (dark) points correspond to light (heavy) systems. 
The solid curve (taken from Ref. [^lj]) approximates the curve for chemical freeze-out of hadronic 
matter. It corresponds to an average energy per hadron of 1 GeV in a hadron-gas model, while the 
crosshatched band indicates location and uncertainty of the boundary between a gas of hadrons 
and a plasma of gluons and deconfined quarks. 
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TABLES 

TABLE I. Mean inclusive ir° multiplicities (M w o) relative to the average number of participat- 
ing nucleons (B) and mean inclusive r] multiplicities (AO relative to (Af^o) measured for various 
nucleus-nucleus collisions in the beam-energy range from 0.8A GeV to 2.0A GeV. The last column 
denotes interpolated mean deuteron multiplicities (M<j) relative to the mean nucleon multiplicity 



(Mff). The data for Kr + Zr and Au + Au (both at 1A GeV, see [18|) have been revised to account 
for the fact that the first-level trigger condition in effect in these measurements introduced a bias 
towards centrality. 



J^beam l^ 1 V J 


UJ o ICL11 


(B\ 
\ D I 


< M ^0> r %1 

(b) 




(M N > L /0 J 


u.o 


V_y ^ K-J 


a 
u 


q 7 _i_ n q 
O. / ZC u.o 


0.31 ±0.11 


ii q _i_ o 7 


0.8 


Ar + Ca 


20 


3.1 ±0.5 


0.41 ±0.04 


16.3 ±3.9 


0.8 


Au + Au 


125±15 


1.6 ±0.3 


0.38 ±0.08 


18.4 ±2.5 


1.0 


C + C 


6 


5.6 ±0.4 


0.57 ±0.14 


10.6 ±2.5 


1.0 


Ar + Ca 


20 


3.0 ±0.3 


1.3 


±0.8 


15.3 ±3.7 


1.0 


Kr + Zr 


79±9 


2.5 ±0.7 


1.3 


±0.6 


17.7 ±2.3 


1.0 


Au + Au 


164±20 


2.3 ±0.5 


1.4 


±0.6 


17.2 ±2.3 


1.5 


Ar + Ca 


20 


6.5 ±0.5 


2.2 


±0.4 


12.7 ±3.0 


2.0 


C + C 


6 


13.8 ±1.4 


3.6 


±0.4 


7.0 ±1.7 


2.0 


Ca + Ca 


20 


11.1 ± 1.1 


2.7 


±0.3 


10.2 ±2.5 


1.9 


Ni + Ni 


29 


8.6 ±0.9 


3.3 


±0.4 


12.7 ±2.8 
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TABLE II. Baryon resonances considered as constituents of the hadronic fireball. Given for 
each resonance are the nominal mass m#, the total width Tr at m# as well as the partial widths 



correspond 


ing to the three decay modes considered 


in the present model. 


Values are 


taken from 


H (3- or 


4-star status) with the exception of A(1232) p3| and 7V(1535) 


& 




resonance 


mass m/j [MeV] 


widtn 1 r [MeVJ 


T /T ro/i T 
I Itt/1 R 1%J I 2tt 


/r« [%] 


r v /r R [%) 


A(1232) 


1232 


110 


100 








iV(1440) 


i A a n 
1441) 


350 


r - r 
DO 


35 





iV(1520) 


152U 


1 on 

120 


55 


45 





iV(1535) 


1 CT A A 

1544 


203 


50 





50 


A(1600) 


1600 


350 


1 r 

15 


85 





A(1620) 


1620 


150 


30 


70 





iV(1650) 


1650 


150 


80 


20 





JV(1675) 


1675 


150 


45 


55 





iV(1680) 


1680 


130 


65 


35 





iV(1700) 


1700 


100 


10 


90 





A(1700) 


1700 


300 


15 


85 





iV(1710) 


1710 


100 


15 


85 





iV(1720) 


1720 


150 


20 


80 
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TABLE III. Particle ratios and chemical freeze-out parameters for peripheral, semicentral, and 
central collisions in the system Au + Au at 0.8A GeV. The impact parameter was selected according 
to the number of participating nucleons (B) given in the first row. The next entries denote the 
meson ratios ((M n o) / (B) , (M r) )/(M 7r o)), the inclusive deuteron/nucleon ratio (Ma)/(M^), and the 
freeze-out radius R c . The subsequent rows show the freeze-out parameters ps and T c together 
with the resulting baryon density ps relative to the nuclear ground-state density po = 0.168 fm~ 3 . 



The values are determined by \ 2 minimization, the uncertainties represent la standard deviations 
as evaluated from the size of the error ellipse of the pb,T c pair at x 2 = Xmin ± 1- 





peripheral 


semicentral 


central 


(B) 


40 ±10 


227 ±20 


345 ± 25 


(M n o)/(B) [%] 


1.3 ±0.4 


1.6 ±0.2 


2.1 ±0.3 


{M v )/{M„o) [%] 


0.27±0.13 


0.29 ±0.10 


0.60 ±0.22 


(M d )/(M N ) [%] 




18.4 ± 2.5 (inclusive) 




Rc [fm ] 


6.5 


11.0 


12.0 


PB [MeV] 


812 ±5 


812 ±5 


803 ±6 


T c [MeV] 


51±2 


52 ±2 


56 ±3 


Pb/po 


0.21 ±0.07 


0.23 ±0.08 


0.28 ±0.09 


Xmin 


0.88 


1.50 


2.93 
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TABLE IV. Chemical freeze-out parameters hb and T c for the various systems investigated, 
together with the resulting baryon densities ps- The values are determined by x 2 minimization, 
uncertainties represent la standard deviations as evaluated from the size of the error ellipse of the 
Pb,T c pair at x 2 = Xmin + 1- Also given are the model-dependent thermal freeze-out temperatures 
Tb and T$r as obtained from Boltzmann fits and within the blast model of Siemens and Rasmussen, 
respectively. 







chemical freeze-out 






thermal freeze-out 


Ebeam 


System 


R c 


PB 


T c 


Pb/po 


X 2 ■ 


T B 


TsR 


[A GeV] 




[fm] 


[MeV] 


[MeV] 






[MeV] 


[MeV] 


0.8 


C + C 


4.5 


758 ±5 


56 ±2 


0.09 ±0.03 


0.36 


45 ±4 


36 ±4 


n 8 


ill t V < 1 


5.5 


780 ±7 


59 ±1 


n i q + n D5 


10 


54 ±2 


43 ±3 


0.8 


Au + Au 


8.5 


one _i_ k 
808 ± 5 


r A i o 

04 ± 2 


27± 08 


5.90 


5 / ± 4 


A r\ i o 

4y ± D 


1.0 


C + C 


4.5 


737 ±5 


62 ±2 


0.10 ±0.03 


1.05 


52 ±3 


48 ±4 


1.0 


Ar + Ca 


5.0 


779 ±7 


60 ±3 


0.21 ±0.08 


1.14 


68 ±3 


52 ±6 


1.0 


Kr + Zr 


7.5 


790 ±7 


60 ±3 


0.27±0.09 


2.30 


70 ±3 


58 ±6 


1.0 


Au + Au 


9.5 


792 ±7 


58 ±4 


0.27 ±0.09 


3.95 


74 ±4 


61 ±6 


1.5 


Ar + Ca 


4.5 


733 ±7 


76±4 


0.30 ±0.10 


4.46 


78 ±2 


65 ±5 


2.0 


C + C 


3.5 


651 ±8 


89±4 


0.21 ±0.07 


2.39 


84 ±2 


71 ±5 


2.0 


Ca + Ca 


4.5 


685 ±9 


86 ±3 


0.30 ±0.08 


0.96 


86 ±2 


72 ±5 


1.9 


Ni + Ni 


4.5 


707 ±9 


87±4 


0.43 ±0.16 


8.06 


88 ±4 


75 ±6 
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